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Abstract 

^ ■ Two-photon exchange contributions to elastic electron-scattering are reviewed. 



The apparent discrepancy in the extraction of elastic nucleon form factors between 
unpolarized Rosenbluth and polarization transfer experiments is discussed, as well 
as the understanding of this puzzle in terms of two-photon exchange corrections. 
Calculations of such corrections both within partonic and hadronic frameworks are 
reviewed. In view of recent spin-dependent electron scattering data, the relation of 
the two-photon exchange process to the hyperfine splitting in hydrogen is critically 
examined. The imaginary part of the two-photon exchange amplitude as can be ac- 
cessed from the beam normal spin asymmetry in elastic electron-nucleon scattering 
is reviewed. Further extensions and open issues in this field are outlined. 



I like a thing simple but it must be simple through complication. Ev- 
erything must come into your scheme, otherwise you cannot achieve real sim- 
plicity. 

Gertrude Stein, "Afterword," What Are Masterpieces 
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1 Introduction 



Elastic electron-nucleon scattering in the one-photon exchange approximation is a time- 
honored tool to access information on the structure of hadrons. Experiments with increas- 
ing precision have become possible in recent years, mainly triggered by new techniques 
to perform polarization experiments at electron scattering facilities. This opened a new 
frontier in the measurement of hadron structure quantities, such as its electroweak form 
factors, parity violating effects, nucleon polarizabilities, transition form factors, or the 
measurement of spin dependent structure functions, to name a few. For example, ex- 
periments using polarized electron beams and measuring the ratio of the recoil nucleon 
in-plane polarization components have profoundly extended our understanding of the 
nucleon electromagnetic form factors (FFs); for recent reviews on nucleon FFs see e.g. 
Refs. [U E], [3]. For the proton, such polarization experiments access the ratio Ge/Gm 
of the proton's electric ( Ge ) to magnetic ( Gm ) FFs directly from the ratio of the 
"sideways" and "longitudinal" polarizations in elastic electron-nucleon scattering as [I], 



P s I 2e G E (Q* 



Here, 



P z Vr(l + e) Gm(Q 2 



t = ^— , - = 1 + 2(1 + r)tan 2 - 



4M 2 
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Q 2 = —q 2 is the momentum transfer squared, 9 is the laboratory scattering angle, and 
< e < 1. Recently, this ratio has been measured at the Jefferson Laboratory out to a 
space-like momentum transfer Q 2 of 5.6 GeV 2 El [7] . It came as a surprise that these 
experiments extracted a ratio of Ge/Gm which is clearly at variance with unpolarized 
measurements [HJ [9j [10] using the Rosenbluth separation technique, which measures the 
angular dependence of the differential cross section for elastic electron-nucleon scattering 
at fixed Q 2 , 

oc G 2 M + £ -G\ , (3) 



dflLab T 

where the proportionality factor is well known, and isolates the e-dependent term. In each 
case, the quoted formulas assume single-photon exchange between electron and nucleon. 

To explain the discrepancy between the two experimental techniques, suspicion falls 
on the Rosenbluth measurements, but not because of experimental problems perse. The 
Rosenbluth formula, Eq. ([3]), at high Q 2 has a numerically big term, G 2 M , and a small 
term. The results for Ge come from the small term. Any omitted ^-dependent corrections 
to the large term can thus have a strikingly large effect on Ge- 

Two-photon exchange, Fig.[H is one thinkable culprit. The subject has a long history. 
In fact as early as the late 1950s and during the 1960s, when electron-nucleon scattering 
was measured systematically at the Stanford Linear Accelerator Center to access nu- 
cleon electromagnetic FFs, the validity of the one-photon exchange approximation has 
been discussed both theoretically and experimentally. Because the one-photon exchange 
cross section depends quadratically on the lepton charge, the difference between electron- 
nucleon and positron-nucleon cross sections is a test for two- or multi-photon exchange 
processes. Early comparisons of electron- and positron-nucleon scattering cross sections 
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were consistent with equal cross sections [HI [12], but the precision achieved in those early 
investigations could not exclude two-photon exchange effects at the few percent level of 
the cross section. Contributions of this size can be expected due to the additional elec- 
tromagnetic coupling in the two-photon exchange diagram, which brings in a suppression 
factor a = e 2 /{An) ~ 1/137. 

On the theory side, the calculation of corrections to elastic electron-nucleon scattering 
of order e 2 relative to the Born approximation, also known as radiative corrections, have 
a long history, see e.g. Refs. [T3[ H4"] . which were heavily applied in the analysis of early 
electron-nucleon scattering experiments. The infrared divergences associated with one 
photon in Fig. [1] being soft, i.e., having a vanishingly small four-momentum, had been 
extracted — they cancel against infrared divergences from soft bremsstrahlung — but the 
contributions when both photons are hard {i.e., when the momentum transfer of both 
individual photons is large) were not calculated in those works because of insufficient 
knowledge of the intermediate hadron state. The early calculators were well aware of the 
omission, and explicitly expressed hope that the missing contributions would be small. 
Some early estimates of the two-photon exchange contribution of Fig. [1] with two hard 
photons were attempted though in the late 1950s by Drell and collaborators [T5l fl6] . 
Those works constructed a non-relativistic model for the blob in Fig. [U including the 
nucleon and lowest nucleon resonance contribution, the A(1232). The calculation found 
that the resonance contribution to the two-photon exchange diagram affects the cross 
sections at the ~ 1 % level. Due to the non-relativistic nature of the calculation, the 
result was limited to about 1 GeV electron beam energy. In subsequent works, e.g. 
Refs. [HI EE], two-photon exchange effects were approximately calculated to higher ener- 
gies. In particular, Greenhut [TS] evaluated the contribution of higher nucleon resonance 
intermediate states, with masses up to 1.7 GeV, when evaluating the blob in Fig. [TJ It 
was found that the dispersive (real) part of the two-photon amplitude yields an electron 
to positron cross section ratio which deviates from unity at the 1-2 % level in the few 
GeV region. The relative smallness of the resonance contribution partly originates be- 
cause the real parts of the resonance amplitudes change sign in the integrand entering 
the evaluation of the box diagram. 

Triggered by the experimental discrepancy between polarization transfer and Rosen- 
bluth measurements of the proton form factor ratio Ge/Gm a t larger Q 2 in recent years, 
the field has seen a new life. In 2003, it was noticed in Ref. [H] that the general form of the 
two-photon exchange graphs could be expressed in an effective current x current form, 
but with an extra structure beyond those that gave Gm and Ge- Further, if this extra 
term had the size one might estimate from perturbation theory, then its interference with 
the one-photon exchange amplitude could be comparable in size to the (Ge) 2 term in the 



Figure 1: Two photon exchange in the lepton - nucleon scattering process l(k) + N(p) 
l(k') + N(p'), with k, k' (p,p') the four-momenta of leptons (nucleons) respectively. 
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Rosenbluth cross section at large Q 2 . In addition, there could be e-dependent modifica- 
tions to the (Gm) 2 term. Hence there was motivation for a precise evaluation. Realistic 
calculations of elastic electron-nucleon scattering beyond the Born approximation are 
required in order to demonstrate in a quantitative way that 27 exchange effects are in- 
deed able to resolve this discrepancy at larger Q 2 . In the current work, we will review 
several such attempts and describe the present status of this field. On the experimental 
side, recent years have seen first attempts to extract the effect of two-photon exchange 
contributions in a quantitative way from the electron-proton scattering data [20] . 

Besides offering a way to explain the glaring discrepancy between two methods of 
measuring the proton electric form factor, the study of two-photon processes also contain 
opportunities to access nucleon structure physics which surpasses the information con- 
tained in nucleon FFs. The possibility arises because a successful two-photon calculation 
involving a hadronic system requires knowledge of hadronic structure, of a sort which 
has only been available recently. For example, one line of investigation arises when the 
virtuality of one or both of the photons in the two-photon process is large compared to 
a nucleon mass scale. In that case, the hard scale allows one to access the Compton 
scattering subprocess on a quark within the nucleon. The new (non-perturbative) pieces 
of information which one then accesses from such a process are the quark correlation 
functions within the nucleon, also known as generalized parton distributions (GPDs); for 
recent reviews see Refs. [2T| |22| [23l [2H [25]. We will review how these GPDs can in turn 
be used to estimate the two-photon exchange diagram of Fig. [I] at large Q 2 . 

Another line of calculations involves the hadronic corrections to ultra-high precision 
atomic physics experiments, such as the hydrogen hyperfine-splitting. Theoretical pre- 
dictions for quantities such as the Lamb shift in hydrogen or the hydrogen hyperfine 
splitting can currently be performed within quantum electrodynamics to such accuracy 
that the leading theoretical uncertainties are related to the nuclear size or the nuclear 
excitation spectrum; for reviews see e.g. Refs. [261 [27]. For the hyperfine splitting in hy- 
drogen, which is known to 13 significant figures, current theoretical understanding is at 
the part per million level. The leading theoretical uncertainty involves the calculation of 
the two-photon exchange graph of Fig. [1] for zero momentum transfer between the bound 
electron and the proton, allowing for all nucleon excited states in the blob. The current 
theoretical understanding of these two-photon exchange corrections will be reviewed in 
this work. 

There are several physical problems where a one-photon exchange potential is not suf- 
ficiently accurate. Besides the description of simple atoms, such as hydrogen and helium, 
also for a precise description of positronium one needs to include two- and multi-photon 
exchange effects. In particular in the interaction of electrically neutral systems, such as 
neutral atoms and molecules, the effect of two-photon exchange gives the dominant con- 
tribution to the forces between such systems, for a theoretical review of such dispersion 
forces; see e.g. Ref. [28] . 

To push the precision frontier further in electron scattering as well as in the hadronic 
corrections to atomic physics quantities, one needs a good control of 27 exchange mech- 
anisms and needs to understand how they may or may not affect different observables. 
This justifies a systematic study of such 27 exchange effects, both theoretically and ex- 
perimentally. Besides the real (dispersive) part of the 27 exchange amplitude, which can 
be accessed by reversing the sign of the lepton charge, also precise measurements of the 
imaginary part of the 27 exchange amplitude became possible in very recent years. The 
imaginary (absorptive) part of the 27 exchange amplitude can be accessed through a sin- 
gle spin asymmetry (SSA) in elastic electron-nucleon scattering, when either the target 
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or beam spin are polarized normal to the scattering plane, as has been discussed some 
time ago [291130] . As time reversal invariance forces this SSA to vanish for one-photon ex- 
change, it is of order a = e 2 /(47r) ~ 1/137. Furthermore, to polarize an ultra-relativistic 
particle in the direction normal to its momentum involves a suppression factor m/E (with 
m the mass and E the energy of the particle), which typically is of order 10 -4 — 1CT 3 
for the electron when the electron beam energy is in the 1 GeV range. Therefore, the 
resulting target normal SSA can be expected to be of order 10~ 2 , whereas the beam 
normal SSA is of order 10~ 6 — 10 -5 . Measurements of small asymmetries of the order 
ppm are quite demanding experimentally, but have been performed in very recent years, 
and will also be reviewed in this work. 

The outline of the present work is as follows. In Section 2, we review the elastic 
electron-nucleon scattering beyond the Born approximation and highlight the discrep- 
ancy in the extraction of Ge/Gm using polarization transfer and unpolarized (Rosen- 
bluth) measurements. We give a brief review of the different attempts which have been 
made recently to explain this difference in terms of two-photon exchange corrections, and 
present in more detail a partonic description at larger momentum transfers. 

In Section 3, we describe the hadronic corrections of the hydrogen hyperfme splitting, 
based on the latest evaluation of the forward polarized structure functions which enter 
in the calculation of the two-photon exchange diagram. 

In Section 4, we review the beam and target single spin asymmetries which measure 
the imaginary part of the two-photon exchange amplitude. In particular, we give an 
overview of the recent high precision measurements in case of a polarized beam with 
normal beam spin polarization. 

We conclude in Section 5, and spell out a few open issues in this field. 
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2 Two-photon in elastic electron-nucleon scattering 



2.1 Hard two- phot on exchange 

In this section, we study attempts to evaluate completely the two-photon exchange con- 
tributions to electron-proton elastic scattering, specifically including the exchange of two 
hard photons, which can probe well inside the proton and which require detailed knowl- 
edge of proton structure to evaluate. The immediate motivation, already given in the 
introduction, is the conflict between the Rosenbluth, with pre-2003 radiative corrections, 
and the polarization measurements of Ge/Gm for the proton. The difference between 
the two techniques is a factor 4 at Q 2 = 5.6 GeV 2 ! (The data will appear on plots later 
in this section, when we discuss how well the proposed resolutions of the conundrum are 
working. ) 

Modern quantitative calculations either treat the hadronic intermediate state in Fig. [1] 
as a proton plus a set of resonances, or else treat it in a constituent picture using gener- 
alized parton distributions. 

In the earliest of the modern calculations, Blunden, Melnitchouk, and Tjon [31] eval- 
uated the two-photon exchange amplitude keeping just the elastic nucleon intermediate 
state. They found that the two-photon exchange correction with an intermediate nu- 
cleon has the proper sign and magnitude to partially resolve the discrepancy between 
the two experimental techniques. Later, the same group, joined by Kondratyuk [33J, 
included contributions of the A (1232) in the intermediate state, which partly canceled 
the elastic terms. Most recently, Kondratyuk and Blunden [M] included five more baryon 
resonances in the intermediate state. While finding that the overall contribution of the 
additional resonances was not large, the totality of their corrections with their choices 
for the 7-nucleon-resonance vertices leads to good agreement with the Rosenbluth data 
using the form factors obtained from polarization data. 

Borisyuk and Kobushkin [35] also considered two-photon corrections with elastic nu- 
cleon intermediate states, and used dispersive techniques to reduce the necessary integrals 
to ones involving the vertex form factors with only spacelike momentum transfers. They 
were able to reduce it further to a single numerical integral for sufficiently low Q 2 . (They 
do not show any Rosenbluth-type plots, but their plotted results for (say) SGm/Gm, in 
notation defined below, are in line with results known from what will be described in the 
rest of this section, despite the rather different methodology.) 

In [311 a g rou P including the present authors calculated the hard two-photon 
elastic electron-nucleon scattering amplitude at large momentum transfers by relating 
the required virtual Compton process on the nucleon to generalized parton distribu- 
tions (GPDs), which also enter in other wide angle scattering processes. This approach 
effectively sums all possible excitations of inelastic nucleon intermediate states. It was 
found that the two-photon corrections to the Rosenbluth process indeed can substantially 
reconcile the two ways of measuring Ge/Gm- 

Rosenbluth data is also available where the recoiling proton, rather than the electron, 
is detected [10J. These data appear to match the data where the scattered electron 
was detected. The two-photon exchange contributions are the same whatever particle is 
detected. However, the bremsstrahlung corrections are different. We shall defer detailed 
discussion of the proton-detected data pending reassessment of the original [TH [38] and 
the new [39l [TO] proton-observed bremsstrahlung calculations. 
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2.2 Elastic electron-nucleon scattering observables 

In order to describe elastic electron-nucleon scattering, 

l(k, h) + N(p, A at) - l(k', ti) + N(p', X' N ), (4) 

where h, h', A^, and are helicities, we adopt the definitions 

p + p f k + k i 

P= , K = — ^— , q = k-k = p -p , (5) 

define the Mandelstam variables 

s = {p + k)\ t = q 2 = -Q\ u = (p-k') 2 , (6) 

let v = K ■ P, and let M be the nucleon mass. 
The T-matrix helicity amplitudes are given by 

T^ h XN = (k',ti;p',\' N \T\k,h;p,\ N ) . (7) 

Parity invariance reduces the number of independent helicity amplitudes from 16 to 8. 
Time reversal invariance further reduces the number to 6 [3D]. Further still, in a gauge 
theory lepton helicity is conserved to all orders in perturbation theory when the lepton 
mass is zero. We shall neglect the lepton mass. This finally reduces the number of 
independent helicity amplitudes to 3, which one may for example choose as 

T +,+ . T +,+ . T +,+ = T +,+ _ (g) 

(The phase in the last equality is for particle momenta in the xz plane.) 

Alternatively, one can expand in terms of a set of three independent Lorentz struc- 
tures, multiplied by three generalized form factors. One such T-matrix expansion is 

e 2 _ 

T h ,\' N x N = -^uik'^^k^h) (9) 



Q 

x u(p\\' N ){G M ^-F 2 — + F 3 ^ 1 ^)u(p,\ 



N) 



The expansion is general. The overall factors and the notations Gm and F 2 have been 
chosen [19] to have a straightforward connection to the standard form factors in the 
one-photon exchange limit. 

If desired, one may replace the F3 term by an axial-like term using the identity, 

s — u 

n(/c')7 ■ Pu(k) x u(j/)7 • Ku(p) = — - — u{k')^^u{k) x u(p')^u(p) 

+ -u(k') W u(k) x uip'h^uip) , (10) 

which is valid for massless leptons and any nucleon mass. We will, however, continue 
with the T-matrix in the form shown in Eq. (Q. An equivalent expansion has also been 
studied in Ref . HI] • 
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The scalar quantities Gm, F2, and ^3 are complex functions of two variables, say v 
and Q 2 . We also use 

G e = G m -(1 + t)F 2 . (11) 

To separately identify the one- and two-photon exchange contributions, we use the no- 
tation Gm = Gm + 8Gm, and Ge = Ge + SG E , where Gm and Ge are the usual proton 
magnetic and electric form factors, which are functions of Q 2 only and are defined from 
matrix elements of the electromagnetic current. The amplitudes F3 = 5F^, 5Gm, and 
SG E , originate from processes involving the exchange of at least two photons, and are of 
order e 2 (relative to the factor e 2 in Eq. (jHJ)). 
The unpolarized cross section is 

da to r da NS 



dVt Lah e(l + r) dVL Lab ' 
where the "no structure" cross section is 

da NS _ 4a 2 cos 2 (6/2) E' 3 
dVL Lab ~ E 



(13) 



and E and E' are the incoming and outgoing electron Lab energies. Other quantities 
are defined after Eq. (JT]). The reduced cross section including the two-photon exchange 
correction becomes 



or = G 2 M + £ -G 2 E +2 G M K (SG M + + 2^-Ge K (sG e + ^F 3 ) + (9(e 4 ), (14) 

where TZ stands for the real part. 

The general expressions for the double polarization observables, including two-photon 
exchange, are 



A s = (15) 
'2e(l-e) (2h e ) ( n , n _ / ~ x , n _ / ~ , v ~ 



{g e G m + G E n{5G M ) + G M n (SG E + -^h) + ^(e 4 )} 



Pi = -A = ^^ [G 2 M + 2GMn{sG M + T ^- £ ^F^+0(e 4 )], 

where h e = ±1/2 is the helicity of the electron and we assumed m e = 0. The polarizations 
are related to the analyzing powers A s or Ai, by time- reversal invariance. That the 
polarization Pi is unity in the backward direction, e = 0, follows generally from lepton 
helicity conservation and angular momentum conservation. 



2.3 Two-photon exchange at the quark level 

In order to estimate the two-photon exchange contribution to Gm, E 2 and F 3 at large 
momentum transfers, we will consider a partonic calculation illustrated in Fig. [2j To 
begin, we calculate the subprocess on a quark, denoted by the scattering amplitude H 
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Figure 2: Handbag approximation for the elastic lepton-nucleon scattering at large mo- 
mentum transfers. In the partonic scattering process (indicated by H) , the lepton scatters 
from quarks in the nucleon, with momenta p q and p' . The lower blob represents the GPDs 
of the nucleon. 



k k' k k' 




Figure 3: Direct and crossed box diagrams to describe the two-photon exchange contri- 
bution to the lepton-quark scattering process, corresponding with the blob denoted by 
H in Fig. 

in Fig. [2j Subsequently, we shall embed the quarks in the proton as described through 
the nucleon's generalized parton distributions (GPDs). 
Elastic lepton-quark scattering, 

Kk) + q( Pq ) ^ l(k') + q(p' q ) , (16) 

is described by two independent kinematical invariants, s = (k + p q ) 2 and Q 2 = —t = 
— {k—k') 2 . We also introduce the crossing variable u = (k—p' q ) 2 , which satisfies s+u = Q 2 . 
The T-matrix for the two-photon part of the electron-quark scattering can be written as 

H h ,x = { -^u(k',h) 1 ,u(k,h) ■ u(p' q ,X)(fil fl + ~hl-KP^)u{p q ,X), (17) 

with P q = ip q +p')/2, where e q is the fractional quark charge (for a flavor q), and 
where u(p q , A) and u{p' , A) are the quark spinors with quark helicity A = ±1/2, which 
is conserved in the scattering process for massless quarks. Quark helicity conservation 
leads to the absence of any analog of F2 in the general expansion of Eq. ([91). 

The partonic scattering helicity amplitudes H h \ of Eq. fTTTj) at order 0(e 4 ) are given 
by the two-photon exchange direct and crossed box diagrams of Fig. [3j The two-photon 
exchange contribution to the elastic electron-scattering off spin 1/2 Dirac particles was 
first calculated in Ref. [12], which was verified explicitly for the work reported in [361 137] . 
The amplitude fx, but not ^3, has an infrared (IR) divergence, which we isolate into a 
soft part, i.e., f\ = /1 + f^ ard . The soft part corresponds with the situation where 
one of the photons in Fig. [3] carries zero four-momentum, and is obtained by replacing 
the other photon's four-momentum by q in both numerator and denominator of the loop 
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integral 



This yields, 

n (f? ft 



47T 2 



In 



X 2 



V- 



ln 



■su , 



7T 

y 



where A is an infinitesimal photon mass controlling the IR divergence. The remaining fi 
can be found in [361 E7] . 

The full electron-quark elastic cross section, using Eq. (fT4l) . is 



da 



da 



It 



2H[f: 



21 



dcri 7 (1 + ^27) 



(19) 



where da^ is the cross section in the one-photon exchange approximation, e = —2 su / (s 2 + 
u 2 ) in the massless limit, and one can easily obtain <5 27 from the fa. 



2.4 Calculation using generalized parton distributions 



Having calculated the partonic subprocess, we next discuss how to embed the quarks 
in the nucleon. We begin by discussing the soft contributions. The handbag diagrams 
discussed so far have both photons coupled to the same quark. There are also contri- 
butions from processes where the photons interact with different quarks. One can show 
that the IR contributions from these processes, which are proportional to the products of 
the charges of the interacting quarks, added to the soft contributions from the handbag 
diagrams give the same result as the soft contributions calculated with just a nucleon 
intermediate state [Hj. Thus the low energy theorem for Compton scattering is satisfied. 
In the handbag approximation, the hard parts which appear when the photons couple to 
different quarks, the so-called cat's ears diagrams, are neglected. 

For the real parts, the IR divergence arising from the direct and crossed box diagrams, 
at the nucleon level, is cancelled when adding the bremsstrahlung contribution from the 
interference of diagrams where a soft photon is emitted from the electron and from the 
proton. This provides a radiative correction term from the soft part of the boxes plus 
electron-proton bremsstrahlung which added to the lowest order term may be written as 



a 



soft — 



()■ 



27, soft 



rep > 
u brems/ 



(20) 



where 0"i 7 is the one-photon exchange cross section. In Eq. (1201) . the soft-photon contri- 
bution due to the nucleon box diagram is given by 



^27, soft 




A 2 



^J(s-M 2 )\u 



In 2 



M 2 




(21) 



+ n 



In 2 




where L is the Spence function defined by 

Liz) = - 



t 



(22) 



The bremsstrahlung contribution where a soft photon is emitted from an electron and 
proton line (i.e., by cutting one of the (soft) photon lines in Fig. [3]) was calculated in 
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Ref. |45j . which we verified explicitly, and is for the case that the outgoing electron is 
detected, 



u brems 



2ir 2 



In 



M 



2\2* 



X 2 y(u-M 2 y 



In 



M 2 



M 2 



u 




(23) 



where AE = E' e el — E' e is the difference of the measured outgoing electron lab energy 
(E e ) from its elastic value (E e el ), and y 



't + y/1 + t) 2 . One indeed verifies that 
the sum of Eqs. ( !2Tl[23l) is IR finite. When comparing with elastic ep cross section data, 
which are usually radiatively corrected using the procedure of Mo and Tsai, Ref. [13], we 



have to consider only the difference of our 5 27)SO /t + ^ ems relative to the 0{Z ) part, in 
their notation, of the radiative correction in [13]. Except for the 7r 2 /2 term in Eq. (l2~Tj) . 
this difference was found to be below 10~ 3 for all kinematics considered in Fig. |H 

After some algebra, one obtains the hard 27 exchange contributions to SGm, SGe, 
and E 3 as, 



Xffhard 



^Qhard 



1 + e 
2s 



B. 



(24) 
(25) 

(26) 



with 



.4 

B 

C 



1 dx 



(s-u)ff ard -suf 3 ] 



ix (s — u) 

1 dx\{s-u)~ft ard -su~h 



E 4 (h 9 + Eq 



(s-u) 

£^f^ d sgm(x)J24H q . 



1 x 
1 (ix ~. 



(H«-tE«) 



(27) 



The functions i/ 9 , E q , and if 9 are the generalized parton distributions, which describe 
removing a quark of a certain momentum from a hadron, and replacing it with a quark 
of another momentum (or, even, in a more general case, with one of a different flavor). 
One can see from Fig. [2] that these are just what we need to knit the amplitudes for 
electron-quark scattering into an amplitude involving the proton. There is a fourth 
GPD, E q , that does not enter the two-photon exchange expressions. The GPDs can be 
measured in deeply virtual or wide-angle Compton scattering, and have been reviewed 
in [2H U2\ ESI [21] and elsewhere. The GPD models we use here are detailed in [33 |4"6"] ; 
so also [4Tj . 

From the integrals A, B, and C, and the usual form factors, we can directly construct 
the observables. The cross section is 



CR — CTR,soft + CTR,hard 



{21 
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where 



(JHMrd = (l+e)G M n (A) + J2e(l + e)-G E TZ(B) + (l-e) G M K (C) . (29) 

v T 



From Eqs. ( !20j) to ( 123]) and the discussion surrounding them, we learned that to a good 
approximation the result for the soft part can be written as 

(?R,soft = <T R ,i<y (l + vra + 5 MT ) , (30) 

where 5 MT is the correction given in Ref . [13] . Since the data is very commonly corrected 
using [13], let us define (J R MT corr = <Jr/(1 + 5 mt ). Then an accurate relationship between 
the data with Mo-Tsai corrections already included and the form factors is 

<y R MT co "' = (g 2 m + £ -G\) (1 + na) + a RMrd , (31) 

where the left hand side is what experimenters often quote as radiatively corrected data. 
Since the Mo-Tsai corrections are so commonly made in experimental papers before 
reporting the data, the "MT corr" superscript will be understood rather than explicit 
when we show cross section plots below. The extra terms on the right-hand-side come 
from new two-photon exchange corrections. The reader may marginally improve the 
expression by including with the (1 + na) factor the circa 0.1% difference between our 
actual soft results and those of [H]. Finally, before discussing polarization, the fact 
that a 7r 2 /2 term, or (na) term after multiplying in the overall factors, sits in the soft 
corrections has to do with the specific criterion we used, that of Ref. [43] , to separate the 
soft from hard parts. The term cannot be eliminated; with a different criterion, however, 
that term can move into the hard part. 

The double polarization observables of Eqs. ffTol) are given by 



P s = - y 2g(1 r £) ^ \g e G m + G E K{C) + G MS j^£n (B) + C>(e 4 ) } , (32) 
p l = ^T^2 A_ f G 2 + G M n(A + C) + 0(e A )\. (33) 

O-R 1 J 



2.5 Results 
2.5.1 Cross section 

Figure H] shows the reduced differential cross section for electron-proton scattering <tr, 
for two values of Q 2 . There are three items on each graph. One is the data. The next is 
the straight line, which is the result of the I-7 exchange calculation using Ge/Gm taken 
from the polarization data [6], with a reasonable and commonly used choice for Gm [48] . 
The slope is too flat to fit the data, reflecting the conflict between the polarization 
measurements and the Rosenbluth measurements with the hard 2-7 corrections. Third 
are the slightly curved lines, showing the results of the 2-7 corrections while still using 
the Ge/Gm ratio from the polarization data. Results are shown for two different model 
GPDs, described in [36 1 137 1 146]; they do not greatly differ. (The renormalization of Gm 
that we have allowed does not affect the slope.) One sees that the hard 2-7 corrections 
steepen the average slope and improve the agreement with the data. It is also important 
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Figure 4: Rosenbluth plots for elastic ep scattering: ctr divided by (/ipG^) 2 , with 
Gjy = (1 + Q 2 /0.71)~ 2 . Dotted curves: Born approximation using Ge/Gm from po- 
larization data [3 [6]. Solid curves: full calculation using the modified Regge GPD, for 
the kinematical range — u > M 2 . Dashed curves: same as solid curves but using the 
gaussian GPD. The data are from Ref. [Sj; the Figure and calculation are from Ref. |37j . 

to note the non-linearity in the Rosenbluth plot, which can be checked with a more 
precise experiment. 

Fig.Elpresents the 2-7 results in a different way. The plot shows the extracted Ge/Gm 
vs. Q 2 . One set of data points, falling linearly with Q 2 , is from the polarization exper- 
iments. Another set of data points, roughly constant in Q 2 and plotted with inverted 
triangles, is from Rosenbluth data analyzed using only the Mo-Tsai [13] radiative correc- 
tions. The solid squares show the best fit Ge/Gm from the Ref. [8] data when analyzed 
including the hard 2-7 corrections. One sees that for Q 2 in the 2-3 GeV 2 range, the 
Ge/Gm extracted using the Rosenbluth method including the 2-7 corrections agree well 
with the polarization transfer results; At higher Q 2 , there is at least partial reconciliation 
between the two methods. 

There is one more point on Fig. [51 which shows the result of also including some hard 
bremsstrahlung corrections, which will be discussed below. 

2.5.2 Polarization transfers 

The 2-7 corrections do not impact the polarizations measurements as strongly as the 
Rosenbluth measurements. The left panel of Fig. [6] shows the correction to the P s /Pi 
ratio from the hard 2-7 exchange. Most of the effect is on P s , shown separately in the 
right panel of the same Figure; the effect on Pi is too small to show on Figures like 
these. The present polarization experiments have e ~ 0.7. The 2-7 corrections induce an 
e-dependence that could be seen in a precise experiment [19] . 

2.5.3 Positron-proton vs. electron-proton 

Positron-proton and electron-proton scattering have the opposite sign for the two-photon 
corrections relative to the one-photon terms. Hence one expects e + p and e~p elastic 
scattering to differ by a few percent. Figure [7] shows our results for three different Q 2 
values. These curves are obtained by adding our two-photon box calculation, minus the 
corresponding part of the soft only calculation in [13], to the one-photon calculations; 
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Figure 5: Determinations of the proton Ge/Gm ratio. The polarization data is from 
Gayou et al. [6] and Punjabi et al. [7], and the Rosenbluth data is from Andivahis et al. [8], 
which include only the well-known Mo-Tsai corrections. Our Rosenbluth Ge/Gm include 
the two-photon corrections, and for one point also a hard bremsstrahlung correction, still 
using Andivahis et al. data. Some of our points for the Rosenbluth results are slightly 
offset horizontally for clarity. 



hence, they are meant to be compared to data where the corrections given in [13] have 
already been made. Each curve is based on the gaussian GPD and is cut off at low e when 
—u = M 2 . Early data from SLAC are available [TT]; more precise data are anticipated 
from JLab [50]. (Ref. [UJ used the Meister-Yennie [H] soft corrections rather than those 
of Mo and Tsai. We have checked that for these kinematics the difference between them 
is smaller than 0.1%, which is negligible compared to the size of the error bars.) 



2.5.4 Results from single-baryon intermediate states 

We have focused on a partonic view of the two-photon physics. The results when viewing 
the hadronic intermediate state as a proton, or a proton plus a set of resonances, are 
similar [31j [32j (33[ EH [35]. The effect, in a calculation with just a proton intermediate 
state, of the extra 2-7 corrections upon extracting Ge/Gm from a Rosenbluth experiment 
is shown in Fig. [HI Further, corrections to the polarization experiments are just a few 
percent, in the same direction as found in the partonic evaluation, with nearly all the 
effect coming upon P s and not on Pi [32] . 



2.6 Remarks on related topics 

There has been a suggestion that hard bremsstrahlung may cause the difference between 
the Rosenbluth and polarization results [ST]. Bremsstrahlung means a process where a 
real photon is emitted. If the photon energy is sufficiently low, the experimenters will 
fail to see it and will count the reaction as elastic. Usual bremsstrahlung calculations 
are for soft bremsstrahlung, where the emitted photon energy is kept only to linear order 
in denominators and entirely omitted in numerators. Soft bremsstrahlung multiplies all 
amplitudes by the same factor and does not, for a relevant example, change the slope 
on a Rosenbluth plot. If one makes no approximations in the photon energy, there 
can be different effects on different spin amplitudes. Thus the claim is that emitted 
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2-Y corrections to polarization ratio - proton Sideways polarization - proton 




Figure 6: Polarization P s /Pi and P s compared to the I7 exchange results for ep scattering 
at Q 2 = 5 GeV 2 . The solid and dashed curves show the 2-7 exchange correction using 
the GPD calculation, for the kinematical range where both s, —u > M 2 . The dotted 
curves on the right shows the Born approximation (I7 exchange) result. 
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Figure 7: Ratio of e + /e~ elastic cross sections on the proton. The GPD calculations of 
the 27 correction are for Q 2 of 2, 5, and 9 GeV 2 , for the kinematical range with — u above 
M 2 . Also shown are all known data, from (TTJ, with Q 2 above 1.5 GeV 2 (the missing 
central value is at 1.111). The numbers near the data give Q 2 for that point in GeV 2 . 



photons that are energetic enough to affect the spin structure of the calculation but still 
small enough to escape detection, give rise to the difference between the two methods of 
measuring Ge/Gm- A contrasting numerical claim is that hard bremsstrahlung effects 
are noticeable and helpful in reconciling the Rosenbluth and polarization experiments, 
but are not decisive. Along these lines is a result for hard bremsstrahlung at Q 2 = 6 
GeV 2 from Afanasev [52], which has been added to the 2-7 results and included in Fig. [5j 

These contrasting claims clearly need adjudication, but an independent reexamination 
is not available as of this writing. 

Electron, or muon, scattering off deuterons or larger nuclei has not been within the 
scope of the present review. Larger nuclei have a factor Z advantage in the relative 
size of the 27 and I7 effects, although breakup effects vitiate this advantage for elastic 
scattering except at low energy. One can examine some of the work seeking evidence of 
27 effects in larger nuclei in [531 EH ESI EH EI] • 

Two-photon exchange effects also affect parity- violating e-p elastic scattering via their 
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Figure 8: Extracting Ge/Gm with 27 corrections calculated using a single proton as the 
hadronic intermediate state. "PT" is Ge/Gm obtained from the polarization transfer 
experiments; "LT" is Ge/Gm obtained from a Rosenbluth experiment using only the Mo- 
Tsai radiative corrections, and "LT+27" includes the extra 27 corrections done this way. 
This figure is a based on a figure in Ref. [32], and we thank the authors for providing it. 

interference with the lowest order Z-exchange diagram. Ref. [59] pointed this out, and 
found that the 27 exchange also led to extra terms with different r and e dependences 
than those known from analyses using only Born diagrams. The calculated size of the 
effects, using the partonic model at Q 2 of several GeV 2 , was of (9(1%). This is below 
present experimental uncertainties, but parity-violating experiments with 0((l/2)%) un- 
certainties are planned. 

Arrington and Sick [SB], considering the effects that the most recent and precise 
low-Q 2 determinations of Ge and Gm would have upon parity- violating e-p elastic scat- 
tering, and included the two photon correction terms that were pointed out in [59]. The 
actual two-photon calculations at their Q 2 were done using single hadron intermediate 
states [60] . 
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3 Polarizability corrections to hydrogen hyperfine 
splitting 



3.1 Two-photon exchange and atomic structure 

We begin with some explanation of how this piece of atomic physics fits properly into a 
review of two-photon physics. Hydrogen hyperfine splitting (hfs) in the ground state is 
known to 13 significant figures in frequency units |61j . 

£hfs(e~>) = 1 420.405 751 766 7(9) MHz . (34) 

This accuracy is remarkable to theorists, who are currently hopeful of obtaining a cal- 
culation accurate to a part per million (ppm). To reach this goal, some improvement is 
needed, and the current best calculations are a few ppm away from the data. 

The main uncertainty in calculating the hfs in hydrogen comes from the hadronic, 
or proton structure, corrections. The generic process that contributes is two-photon 
exchange, shown in Fig. [9^, which involves the proton structure because each of the 
exchanged photons could individually be quite energetic. 

One-photon exchange, Fig. [9b, does not involve proton structure, at the accuracy 
needed for the present purpose. The characteristic momentum of the electron in a hy- 
drogen atom is of 0(am e ), which is very low on a nuclear physics scale. Hence the q 2 of 
an exchanged single proton is very low, and the variation of the proton form factor from 
its q 2 = value is minimal. One can show that keeping the electron momentum gives 
corrections of 0(am e /M) smaller than what comes from two-photon exchange. Hence 
one sets the momenta of the electrons to zero. (For information, in the one-photon ex- 
change hfs calculation there comes a q 2 factor in the numerator which cancels the 1/q 2 
from the photon propagator; then the neglect of the electron momenta is safe.) 



e - er 




(a) (b) 



Figure 9: (a) Generic two-photon exchange diagram, giving proton-structure corrections 
to hyperfine splitting, (b) One-photon exchange. 



3.2 Hyperfine splitting calculations 

The calculated hyperfine splitting in hydrogen is [EEJ E21 E3], 

^hfs(e-p) = (1 + A QED + A^ cak + A str ) El ■ (35) 
the two-photon exchange lies in the structure dependent term A str . The Fermi energy is 

%<x>m z r 16a 2 n p 

Ep = nZ. f^B^p = ^- — 7^ 77773 > ( 36 ) 

on 3 [Ib (1 + m e /M) 
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where m r = m e M/(M + m e ) is the reduced mass. By convention, in E P F one uses the 
actual magnetic moment for the proton and the Bohr magneton for the electron (note 
that fiB can be used to replace the electron mass), and Roo is the Rydberg constant 
in frequency units. The second form allows optimal accuracy in evaluating E F . The 
least accurately known quantity is the ratio h p /hb, which is known to 8 figures. Hence 
to the ppm level, E F is known more than sufficiently well. The QED [26l [27] and 
weak interaction corrections [M] are well known and will not be discussed, except to 
mention that the QED corrections could be obtained, for the present purposes, without 
calculation. They are the same as for muonium, so it is possible to obtain them to an 
accuracy more than adequate for the present purpose using muonium hfs data and a 
judicious scaling [65] [66]. 

The structure dependent corrections are, in a standard treatment, divided into Zemach, 
recoil, and polarizability corrections, 

A str = A Z + A P R + A pol . (37) 

The first two terms arise when the blob in Fig. [9^ is just a proton, as in Fig. [TUl 
and they are together called the elastic corrections. The electron-photon vertex is well 
known, and the proton-photon vertex is given by [67J [68] 

= 7^1 (Q 2 ) + ^,uq u F 2 (Q 2 ) (38) 

(for the photon with incoming momentum q) if the intermediate proton is on-shell. Of 
course, it is generally not. However, one can show that the imaginary part of the diagram 
does come only from kinematics where the intermediate electron and proton are on-shell. 
Hence, one can correctly use the above vertex to calculate the imaginary part of the 
diagrams, and then obtain the whole of the diagram using a dispersion relation. 



k k k k 




Figure 10: Two-photon exchange diagrams for the "elastic" proton-structure corrections 
to hyperfine splitting. 



In the nonrelativistic limit, the recoil terms are zero and the Zemach term is not. 
(The nonrelativistic limit is M — > oo, with m e and the proton size held fixed. Proton 
size information is embedded in the form factors F\ and F2.) The Zemach term [69J was 
calculated long ago and in modern form is, 



A: 



8am r [°° dQ 



71 



G E {Q 



1 



K r 



-2am r r z , 



(39) 



the last equality defining the Zemach radius r Z - The charge and magnetic form factors 
are linear combinations of F\ and F2, 



Gm 
Ge 



Fi + F 2 , 

Q 2 



F, 



AM 2 ' 



(40) 
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The recoil corrections won't be explicitly displayed in this paper because they are 
somewhat long (although not awful; see [671 EH])- An important point is that they 
do depend on the form factors and hence upon the proton structure. However, their 
numerical value is fairly steady by present standards of accuracy when they are evaluated 
using different up-to-date analytic form factors based on fits to the scattering data. 

One gets the polarizability corrections by summing all contributions where the inter- 
mediate hadronic state, the blob in Fig. [9^, is not a single photon. Paralleling the elastic 
case, one can show that the imaginary part of this diagram comes only from configura- 
tions where the intermediate electron plus hadronic state is kinematically on-shell, i.e., 
physically realizable. Hence one can calculate the imaginary part of the box diagram if 
one has data on inelastic electron-proton scattering. Then, further paralleling the elastic 
case, one obtains the full box diagram via a dispersion relation. 

To see some detail, the lower half of Fig. [Hh is the same as forward Compton scattering 
of off-shell photons from protons, which is given in terms of the matrix element 

T»„(q,p,S) = ^mJ^ j q< <pS\T {j M (0,j,(0)} \pS) , (41) 

where is the electromagnetic current and the states are proton states of momentum p 
and spin 4-vector S. The spin dependence is in the antisymmetric part 



T A — 1 a n a 



(HM) + HM)) S^-H 2 (u,q 2 



S-q 



p-q 



(42) 



The two structure functions Hi and Hi depend on q 2 and on the lab frame photon energy 
v, defined by Mv — p- q. 

The optical theorem that relates the imaginary part of the forward Compton am- 
plitude to the inelastic scattering cross section for off-shell photons on protons. The 
relations precisely are 



lmHi(u,q 2 ) = - gi(u,q 2 ) and lmH 2 (u, q 2 ) = — g 2 (v, q 2 ) , (43) 
v v & 



where gi and g 2 are functions appearing in the cross section and are measured [701 I7T1 
[TJ [731 at SLAC, HERMES, JLab, and elsewhere. 

Using the Compton amplitude to evaluate the inelastic part of the two-photon loops 
leads to 



Ap i - — 

Hip 

X 



2am R 

(44) 



inel (1 + k p )tt 3 M 

W + d ^ 2 Ql)Q 2 { {2Q2 + Q2 ° )H ' rl( - iQo ' ~ Q2) ~ ^Q 2 Ql H 2%Qo, -Q 2 )} 

where we have Wick rotated the integral so that Qo = —iu, Q = q, and Q 2 = Qq + Q 2 . 
The dispersion relations which gives Hi is, assuming no subtraction, 

where the integral is only over the inelastic region, and a similar relation holds for H 2 . 
The no-subtraction assumption will be discussed later. 
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Integrating what can be integrated, and neglecting m e inside the integral, yields the 
expression [73 EH [77J US ESJ ED] 

^' = W^M (Al + A2h <46) 

where, with r = v 2 /Q 2 , 



A, = 9 - r d S\F2(Q 2 )+^M r^(r)gMQ 2 )) , (47) 



A 2 = -12M r%p 2 (r)g 2 ^Q 2 ). 
Jo Q z Ju th v l 

The auxiliary functions are 



-3r + 2r 2 + 2(2 - tWt(t + 1) 



/3 2 (r) = l + 2r-2^/r(r + l). (48) 

The integral for Ai actually requires further comment. Only the second terms comes 
from the procedure just outlined; it was historically thought convenient to add the first 
term, and then subtract the same term from the the recoil corrections. This stratagem 
allows the electron mass to be taken to zero in A^ The individual terms in Ai diverge 
(they would not had the electron mass been kept), but the whole is finite because of the 
Gerasimov-Drell-Hearn (GDH) [HU E2] sum rule, 

AM — gi ( v , Q ) = -tf (49) 

coupled with the observation that the auxiliary function (3{r) becomes unity as we ap- 
proach the real photon point. 



3.3 Numerics, especially for A po i 

We start this numerical section with a brief discussion of the polarizability corrections. 
They have some history. Considerations of A po i were begun by Iddings in 1965 [75], 
improved by Drell and Sullivan in 1966[76], and given in present notation by de Rafael in 
1971 [77J. But no sufficient spin-dependent data existed, so it was several decades before 
the formula could be evaluated to a result incompatible with zero. In 2002, Faustov and 
Martynenko became the first to use data to obtain results inconsistent with zero [79|. 
They got 

A pol (F&M 2002) = (1.4 ± 0.6) ppm (50) 

However, only SLAC data was available. None of the SLAC data had Q 2 below 0.30 
GeV 2 ; Ai and A 2 are sensitive to the behavior of the structure functions at low Q 2 . 
Also in 2002 there appeared analytic expressions for g lj2 fit to data by Simula, Osipenko, 
Ricco, and Taiuti [S3], which included JLab as well as SLAC data. Simula et al. did not 
integrate their results to obtain A po i, but had they done so, they would have obtained 
A pol = (0.4 ±0.6) ppm [80]. 
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More recently Faustov and Martynenko, joined by Gorbacheva, [83] have reanalyzed 
their result for A po i, and obtained a somewhat larger value, 

A pol (FGM 2006) = (2.2 ± 0.8) ppm . (51) 

The data underlying this result, however, still does not include the lower Q 2 data from 
JLab that will be noted immediately below. 

Data for g±(v, q 2 ) is now improved thanks to the EG1 experiment at JLab, which had 
its first data run in 2000-2001. A preliminary analysis of this data became available in 
2005 [73]; final data is anticipated "soon". The Q 2 range measured in this experiment 
went down to 0.05 GeV 2 . Using analytic forms checked against the preliminary data, 
Ref. [HQ] has evaluated A po i and obtained 

A pol (NCG 2006) = (1.3 ± 0.3) ppm . (52) 

This is similar to the 2002 Faustov-Martynenko result, but with a claim that the newer 
data allows a smaller uncertainty limit. 

A list of the numerical values of the corrections compared to the experimental value 
of the hfs is given is Table [TJ For the polarizability corrections, we used the value from 
Ref. [80J on the grounds that is was based on the most complete inelastic electron-proton 
scattering data. For the Zemach term, we used the value [86] based on the form factor 
fits of Sick [87], because those fits emphasized the low-Q 2 elastic scattering data that 
dominates the Zemach integral. The values for the recoil terms and weak interaction 
corrections have lower uncertainty limits. We took the former come from [65] and they 
are also discussed in [62] ; the latter may be found in [62], [63] • 

Table 1: Up-to-date corrections to hydrogenic hyperfine structure. The first line with 
numbers gives the "target value" based on the experimental data and the best evaluation 
of the Fermi energy (8 figures) based on known physical constants. The corrections are 
listed next. (The Zemach term includes a 1.53% correction from higher order electronic 
contributions [53], as well as a +0.07 ppm correction from muonic vacuum polarization 
and a +0.01 ppm correction from hadronic vacuum polarization [62].) The total of all 
corrections is 1.68 ± 0.60 ppm short of the experimental value. 



Quantity 


value (ppm) 


uncertainty (ppm) 


(E Ms (e~p)/E p F ) 1 


1 103.48 


0.01 


Aqed 


1 136.19 


0.00 


Az (using Friar & Sick [86J) 


-41.59 


0.46 




5.84 


0.15 


A pol (from [80]) 


1.30 


0.30 


^weak 


0.06 




Total 


1101.80 


0.60 


Deficit 


1.68 


0.60 



Thus the sum of what one may argue are the best calculated corrections falls short of 
the data by about 2 ppm, or about 2.8 standard deviations. Of course, some judgement 
has entered the choice of numbers. Other form factor fits to ostensibly the same data give 
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different values for Az', for example the form factors of Kelly [87] lead to Az = — 40.93±49 
ppm. The reader will surely notice that using the Kelly based value of Az and the latest 
value of Ap i from Faustov, Gorbacheva, and Martynenko leads to excellent agreement 
of the measured hfs. 



3.4 Comments on the derivations of the formulas 

The polarizability corrections depend on theoretical results that are obtained using un- 
subtracted dispersion relations. One would like to know just what this means, and since 
there is at least a small discrepancy between calculation and data, one would like to be 
able to assess the validity of such dispersion relations. 

Dispersion relations [88] involve imagining some particular real variable to be a com- 
plex one and then using the Cauchy integral formula to find the functions of that variable 
at a particular point in terms of an integral around the boundary of some region. This 
is a useful thing to do because, at least in particular cases, we have information from 
other sources about the function on part of the boundary, and legitimate reasons for 
neglecting contributions from the rest of the boundary. In the present case we consider 
the functions Hi(is,q 2 ) and we "disperse" in is, treating q 2 as a constant while we do so. 
Three things are needed to make the dispersion calculation work: 

• The Cauchy formula and knowing the analytic structure (locations of poles and 
cuts) of the desired amplitudes. 

• The optical theorem, to relate the forward Compton Imifj to inelastic scattering 
cross sections. 

• A reason for discarding contributions from some oo contour, if the dispersion rela- 
tion is to be "unsubtracted." 

The first two requirements are not in question. 

One can consider the elastic and inelastic parts of Hi separately, but it is best to 
consider them together and separate the terms are the end. For Hi, one can show that 
it is even in is, so we will let the dispersion variable be is 2 rather than is. The contour 
of integration is shown in Fig. [IT], where one should imagine the outside circle having 
infinite radius. Along the real axis, the isolated pole corresponds to elastic scattering and 
the cut is for inelastic scattering kinematics. The result for H\ at some general point v 2 
begins its existence as 



= + i J + J_ / ngJliS . (53) 

Z/gj — IS Z IX Jcut IS' — V A Zlll J\v'\=oo is — is 

The numerator of the first term is the residue (Res) from the poles in v for the elastic 
part of Hi, as from Fig. [TOl 

We can interject here that an alternative calculation of the elastic contributions can 
be done directly, with no dispersion relations, simply using the photon-proton vertex 
given earlier (Eq. (1381 ). whether or not the intermediate proton is on-shell. We do not 
recommend doing the calculation this way, since the vertex cannot be guaranteed correct 
for off-shell protons, but the result is instructive. For purposes of discussion we quote 
the result for Hi. 

x = _2M ( q 2 Fi{q 2 )G M {g 2 ) F%tf)\ 

1 TT \(q 2 + ie) 2 - 4MV 2 AM 2 )' 1 ' 
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Figure 11: Contour in complex v 2 plane for applying Cauchy identity to H\ or H2. 



One can obtain the residue for the dispersion relation from the first term above (and 
do so correctly, since the result is correct at the pole), and insert it into the dispersion 
relation to reproduce (of course) the same term. Of more interest is that the _F| term is 
constant is v, has no imaginary part, and is therefore absent in the dispersive calculation. 
Thus there is a difference, at least at the moment, between the elastic part of the result 
obtained from the dispersive calculation and what one finds in, for example, [67J or |68j . 

The second term leads to the g± term in the quantity Ai given earlier, after using the 
optical theorem to relate 1m Hi to g\. 

The third term is the integral over the part of the contour which is the infinite radius 
circle. The commonly quoted results for A po i, which appear in this paper, depend on 
dropping this term. The term is zero, if Hi falls to zero at infinite \u\. Assuming this is 
true, however, appears to be a dramatic assumption. For example, if the above Hf were 
correct (and it could be: the oft criticized vertex in Eq. (I38j) is not guaranteed to be right, 
but neither are we aware of a guarantee that it is wrong), the assumption would fail for 
Hf alone. Hence, for the assumption to succeed requires an exact cancelation between 
elastic and inelastic contributions, or a failure of Eq. (|54|) on the big contour. On the 
positive side are several considerations. One is that nearly the same derivation gives the 
GDH sum rule, which is checked experimentally and works, within current experimental 
uncertainty (8%) [SH]. Also, the GDH sum rule has been checked in lowest order and 
next-to-lowest order perturbation theory in QED, where it appears to work [901 I9T] - 
Finally, Regge theory suggests the Compton amplitude does fall to zero with energy [92J, 
as one would like, although Regge theory famously gave wrong high v behavior for spin- 
independent analogs of g\ and g^ [93]. On the other hand, parton model calculations [91] 
have suggested a reason why the Regge theory would fail for the spin-idependent structure 
functions but still be correct for the the spin- dependent ones. Hence there are indications, 
though not decisive proof, supporting the unsubtracted dispersion relation. 

The dispersive derivation finishes by subtracting a term involving F 2 2 from the rela- 
tivists recoil term, so as to obtain exactly the elastic corrections A el = A z + A^ that 
were obtained (say) by Bodwin and Yennie for a calculation of the elastic terms only, 
using Eq. (1381) at the photon-proton vertices and no dispersion theory [67] • After adding 
the same term to the polarizability corrections in Ai, one obtains the commonly quoted 
result for Ai [761 123 [79]. As noted earlier, this reorganization also allows Ai to be finite 
in the m e — > limit. Beyond the historical connection, if one is comfortable with the 
unsubtracted dispersion relation, the use of the dispersion theory gives a more secure 
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result because it uses only the pole part of the photon-proton-proton vertex, so that the 
combined elastic and inelastic result does not depend on the general validity of whatever 
photon-proton-proton vertex one uses. 

3.5 Remarks and prospects 

Thus the calculated hyperfme splitting in atomic hydrogen is an example of two-photon 
physics, and requires proton structure information measured at nuclear and particle 
physics laboratories. Until 2006, the largest uncertainty was in the proton polarizability 
corrections, which are related to data from polarized inelastic electron-proton scattering. 
The numerical value of the polarizability contributions to hydrogen hyperfme structure, 
based on latest proton structure function data is A po i = (1.3 ± 0.3) ppm. This is quite 
similar to the Faustov-Martyenko 2002 result, which we think is remarkable given the 
improvement in the data upon which it is based. Most of the calculated A po i comes from 
integration regions where the photon four-momentum squared is small, Q 2 < 1 GeV 2 . 

There is still a modest discrepancy between what we think are the best hydrogen hfs 
calculations and experiment, on the order of 2 ppm. Where may the problem lie? It could 
be in the use of the unsubtracted dispersion relation; or it could be in the value of the 
Zemach radius, which taken at face value now contributes the largest single uncertainty 
among the hadronic corrections to hfs; or perhaps it is a low Q 2 surprise in g\ or #2- It 
is at any rate not a statistical fluctuation in the hfs data itself. 

An interplay between the fields of atomic and nuclear or particle physics may be 
relevant to sorting out the problem. For one example, the best values of the proton 
charge radius currently come from small corrections accurately measured in atomic Lamb 
shift |87j . The precision of the atomic measurement of the proton charge radius can 
increase markedly if the Lamb shift is measured in muonic hydrogen [95], and data may 
be taken in 2007 at the Paul Scherrer Institute. In the present context, the charge radius 
is noticed by its effect on determinations of the Zemach radius. 

We close this section by noting that the final EG1 data analysis from JLab/CLAS 
should be released soon, and this may shift the value of A po i somewhat. We may also 
note that one can keep the lepton masses so as to calculate muonic hydrogen hyperfme 
splitting, and calculations have already appeared [HUES]- 
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4 Beam and target normal spin asymmetries 



In this section we discuss the imaginary part of the two-photon exchange amplitude. It 
can be accessed in the beam or target normal spin asymmetries in elastic electron-nucleon 
scattering, and measures the non-forward structure functions of the nucleon. After briefly 
reviewing the theoretical formalism, we discuss calculations in the threshold region, in 
the resonance region, in the diffractive region, corresponding with high energy and for- 
ward angles, as well as in the hard scattering region. 

The imaginary (absorptive) part of the 27 exchange amplitude can be accessed through 
a single spin asymmetry (SSA) in elastic electron-nucleon scattering, when either the 
target or beam spin are polarized normal to the scattering plane, as has been discussed 
some time ago [291 l30l [97] . As time reversal invariance forces this SSA to vanish for 
one-photon exchange, it is of order a = e 2 /(47r) ~ 1/137. Furthermore, to polarize an 
ultra-relativistic particle in the direction normal to its momentum involves a suppression 
factor m/E (with m the mass and E the energy of the particle), which for the electron 
is of order 10~ 4 — 10 -3 when the electron beam energy is in the 1 GeV range. There- 
fore, the resulting target normal SSA can be expected to be of order 10~ 2 , whereas the 
beam normal SSA is of order 10 -6 — 10~ 5 . A measurement of such small asymmetries is 
quite demanding experimentally. However, in the case of a polarized lepton beam, asym- 
metries of the order ppm are currently accessible in parity violation (PV) elastic eN 
scattering experiments. The parity violating asymmetry involves a beam spin polarized 
along its momentum. However the SSA for an electron beam spin normal to the scatter- 
ing plane can also be measured using the same experimental set-ups. First measurements 
of this beam normal SSA at beam energies up to 1 GeV have yielded values around —10 
ppm [98l l99l 1100] in the forward angular range and up to an order of magnitude larger 
in the backward angular range |101j . At higher beam energies, first results for the beam 
normal SSA in elastic electron-nucleon scattering experiments have also been reported 
recently [HM IT021, ITUB] . 

First estimates of the target normal SSA in elastic electron-nucleon scattering have 
been performed in [301 [97] . In those works, the 27 exchange with nucleon intermediate 
state (so-called elastic or nucleon pole contribution) has been calculated, and the inelas- 
tic contribution has been estimated in a very forward angle approximation. Estimates 
within this approximation have also been reported for the beam normal SSA [104j . The 
general formalism for elastic electron-nucleon scattering with lepton helicity flip, which is 
needed to describe the beam normal SSA, has been developed in [105] . Furthermore, the 
beam normal SSA has also been estimated at large momentum transfers Q 2 in [105j using 
a parton model, which was found crucial [36] to interpret the results from unpolarized 
electron-nucleon elastic scattering, as discussed in Section 2. In the handbag model of 
Refs. [36ll3Tl[105j . the corresponding 27 exchange amplitude has been expressed in terms 
of generalized parton distributions, and the real and imaginary part of the 27 exchange 
amplitude are related through a dispersion relation. Hence in the partonic regime, a 
direct comparison of the imaginary part with experiment can provide a very valuable 
cross-check on the calculated result for the real part. 

To use the elastic electron-nucleon scattering at low momentum transfer as a high 
precision tool, such as in present day PV experiments, one may also want to quantify 
the 27 exchange amplitude. To this aim, one may envisage a dispersion formalism for 
the elastic electron-nucleon scattering amplitudes, as has been discussed some time ago 
in the literature [56l [55] . To develop this formalism, the necessary first step is a precise 
knowledge of the imaginary part of the two-photon exchange amplitude, which enters 
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Figure 12: The two-photon exchange diagram. The filled blob represents the response 
of the nucleon to the scattering of the virtual photon. In the imaginary part of the 
two-photon amplitude, the intermediate state indicated by the vertical dashed line is 
on-shell. 



in both the beam and target normal SSA. Using unitarity, one can relate the imaginary 
part of the 27 amplitude to the electro-absorption amplitudes on a nucleon, see Fig. [121 



4.1 Single spin asymmetries in elastic electron- nucleon scatter- 
ing 

An observable which is directly proportional to the imaginary part of the two- (or multi-) 
photon exchange is given by the elastic scattering of an unpolarized electron on a proton 
target polarized normal to the scattering plane (or the recoil polarization normal to 
the scattering plane, which is exactly the same assuming time-reversal invariance). For 
a target polarized perpendicular to the scattering plane, the corresponding single spin 
asymmetry, which we refer to as the target normal spin asymmetry (A n ), is defined by : 



4. = ^, (55) 

where cr-f (crj_) denotes the cross section for an unpolarized beam and for a nucleon spin 
parallel (anti-parallel) to the normal polarization vector, defined as : 

S£ =((),£„), S n = (kxk')/\kxk'\. (56) 

As has been shown by de Rujula et al. [SU], the target (or recoil) normal spin asym- 
metry is related to the absorptive part of the elastic eN scattering amplitude as : 

, 2 lm(£ spina TZy ■ Abs T 27 ) 
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where T\ 7 denotes the one-photon exchange amplitude. Since the one-photon exchange 
amplitude is purely real, the leading contribution to A n is of order 0(e 2 ), and is due to an 
interference between one- and two-photon exchange amplitudes. When neglecting terms 
which correspond with electron helicity flip (i.e. setting m e = 0), A n can be expressed in 
terms of the invariants for electron-nucleon elastic scattering, defined in Eq. ([9]) as [36] : 
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where X denotes the imaginary part. 

For a beam polarized perpendicular to the scattering plane, one can also define a 
single spin asymmetry, analogously as in Eq. ( l55l) as noted in Ref. |104j . where now a\ 
(crj) denotes the cross section for an unpolarized target and for an electron beam spin 
parallel (anti-parallel) to the normal polarization vector, given by Eq. ( 1561) . We refer to 
this asymmetry as the beam normal spin asymmetry (B n ). It explicitly vanishes when 
m e = as it involves an electron helicity flip. The general electron- nucleon scattering 
amplitude including lepton helicity flip involves six invariant amplitudes and has been 
worked out in Ref. [105] . where the expression for B n can also be found. As for A n , also 
B n vanishes in the Born approximation, and is therefore of order e 2 . 

4.2 Imaginary (absorptive) part of the two-photon exchange 
amplitude 

In this section we discuss the relation between the imaginary part of the two-photon 
exchange amplitude and the absorptive part of the doubly virtual Compton scattering 
tensor on the nucleon, as shown in Fig. [T2l 

The discontinuity of the two-photon exchange amplitude, shown in Fig. [12], can then 
be expressed as : 

——±—u(k', ti^irh+me^Mk, hy-p^-W^tf, X' N ;p, X N ) , (59) 

where the momenta are defined as indicated on Fig. [12], with q x = k — ki, q 2 = k' — ki, and 
Qi — Q2 = Q- Here h{h') denote the helicities of the initial (final) electrons and Ajv(A^) 
denote the helicities of the initial (final) nucleons. In Eq. (1591) . the hadronic tensor 
W^ u (p', \' N ; p, Ayv) corresponds with the absorptive part of the doubly virtual Compton 
scattering tensor with two space-like photons : 

w^(p',\' N]P ,\ N ) = J2W 64 (p + i^-px) < P ',\' N \J^mx >< x\r(o)\ P ,x N >, 

(60) 

where the sum goes over all possible on-shell intermediate hadronic states X (denoting 
p\ = W 2 ). Note that in the limit p' = p, Eq. ([60]) reduces to the forward tensor for 
inclusive electron-nucleon scattering and can be parametrized by the usual 4 nucleon 
forward structure functions. In the non-forward case however, the absorptive part of the 
doubly virtual Compton scattering tensor of Eq. (I60p which enters in the evaluation of 
target and beam normal spin asymmetries, depends upon 18 invariant amplitudes [106j . 
Though this may seem as a forbiddingly large number of new functions, we may use the 
unitarity relation to express the full non-forward tensor in terms of electroproduction 
amplitudes 7*iV — > X. The number of intermediate states X which one considers in 
the calculation will then put a limit on how high in energy one can reliably calculate 
the hadronic tensor Eq. (1601) . In the following section, the tensor will be discussed 
for the elastic contribution (X = N), in the resonance region as a sum over all ttN 
intermediate states (i.e. X = ttN), using a phenomenological state-of-the-art calculation 
for the 7*iV — > -kN amplitudes, in the diffractive region (corresponding with high energy, 
forward scattering) where it can be related to the total photo-absorption cross section 
on a proton, as well as in the hard scattering region where it can be related to nucleon 
generalized parton distributions. 
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There are special regions in the phase space integral of Eq. f|59l) . corresponding with 
near singularities, which may give important contributions (logarithmic enhancements) 
under some kinematical conditions. When the intermediate and initial electrons are 
collinear, then also the photon with momentum q\ is collinear with this direction. For 
the elastic case (W = M) this precisely corresponds with the situation where the first 
photon is soft (i.e. qi — > 0) and where the second photon carries the full momentum 
transfer Q\ ~ Q 2 . For the inelastic case (W > M) the first photon is hard but becomes 
quasi-real (i.e. Q\ ~ rn 2 e )- In this case, the virtuality of the second photon is smaller than 
Q 2 . An analogous situation occurs when the intermediate electron is collinear with the 
final electron. These kinematical situations with one quasi-real photon and one virtual 
photon correspond with quasi virtual Compton scattering (quasi- VCS). Besides the quasi- 
VCS singularities, the two-photon exchange amplitude also has a near singularity when 
the intermediate electron momentum is soft (i.e. |A;i| — > 0). In this case both photons 
are hard but have virtualities which become very small, and vanish if the electron mass 
is taken to zero. This situation, with two quasi-real photons, occurs when the invariant 
mass of the hadronic state takes on its maximal value W max = y^s — m e , and corresponds 
to quasi-real Compton scattering (quasi- RCS). 

4.3 Results and discussion 

4.3.1 Threshold region 

In Ref. |107j . the beam normal spin asymmetry was studied at low energies in an ef- 
fective theory of electrons, protons and photons. This calculation, in which pions are 
integrated out, effectively corresponds with the nucleon intermediate state contribution 
only, expanded to second order in E e /M. To this order, the calculation includes the 
recoil corrections to the scattering from a point charge, the nucleon charge radius, and 
the nucleon isovector magnetic moment. One sees from Fig. [13] (right panel) that the 
theory expanded up to second order in E e /M (indicated by the full results) is able to 
give a good account of the SAMPLE data point at the low energy E e = 0.2 GeV. 

However, when doing the full calculation for the iV intermediate state, which is model 
independent (as it only involves on-shell ^y*NN matrix elements), the result is further 
reduced as seen in Fig. [13] (left panel). Inclusion of threshold pion electroproduction 
contributions, arising from the irN intermediate states, partly cancels the elastic contri- 
butions. Because in this low-energy region, the matrix elements are rather well known, 
it is not clear at present how to get a better agreement with the rather large asymmetry 
measured by SAMPLE [HE]- 

4.3.2 Resonance region 

When measuring the imaginary part of the elastic eN amplitude through a normal SSA 
at sufficiently low energies, below or around two-pion production threshold, one is in a 
regime where these electroproduction amplitudes are relatively well known using pion 
electroproduction experiments as input. As both photons in the 27 exchange process 
are virtual and integrated over, an observable such as the beam or target normal SSA 
is sensitive to the electroproduction amplitudes on the nucleon for a range of photon 
virtualities. This may provide information on resonance transition form factors comple- 
mentary to the information obtained from current pion electroproduction experiments. 
In Ref. |108j . the imaginary part of the two-photon exchange amplitude was calculated 
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Figure 13: Beam normal spin asymmetry for e~*p — > e~p at a beam energy E e = 0.2 GeV 
as function of the cm. scattering angle. Left panel : calculation from Ref. |108j for 
different hadronic intermediate states (X) in the blob of Fig.[T2l : N (dashed curve), ttN 
(dashed-dotted curve), sum of the N and ttN (solid curve). Right panel : calculation 
of Ref. |107|. where the nucleon intermediate state is expanded to leading order (dotted 
curve) and next to leading order (solid curve) in E e /M. For comparison, also the full 
nucleon intermediate state result (dashed curve, same as on left panel) is shown. The 
data point is from the SAMPLE Coll. [HE]. 



by relating it through unitarity to the contribution of X = N and X = ttN intermedi- 
ate state contributions. For the ttN intermediate state contribution, the corresponding 
pion electroproduction amplitudes were taken from the phenomenological MAID analy- 
sis |109j . which contains both resonant and non-resonant pion production mechanisms. 
The calculation of [108] shows that at forward angles, the quasi-real Compton scattering 
at the endpoint W = W max only yields a very small contribution, which grows larger 
when going to backward angles. This quasi-RCS contribution is of opposite sign as the 
remainder of the integrand, and therefore determines the position of the maximum (ab- 
solute) value of B n when going to backward angles. 

In Fig. dH the results for B n are shown at different beam energies below E e — 1 GeV. 
It is clearly seen that at energies E e = 0.3 GeV and higher the nucleon intermediate 
state (elastic part) yields only a very small relative contribution. Therefore B n is a di- 
rect measure of the inelastic part which gives rise to sizeable large asymmetries, of the 
order of several tens of ppm in the backward angular range, mainly driven by the quasi- 
RCS near singularity. First results from the A4 Coll. for B n at backward angles (for E e 
around 0.3 GeV) indeed point towards a large B n value of order -100 ppm for 9 cm around 
150 deg [101] . At forward angles, the sizes of the predicted asymmetries are compatible 
with the first high precision measurements performed by the A4 Coll. [99J, though the 
model slightly overpredicts (in absolute value) B n at E e = 0.570 GeV and 0.855 GeV. 

4.3.3 High-energy, forward scattering (diffractive) region 

At very high energies and forward scattering angles (so-called diffractive limit), it was 
shown in Refs. [1101 1111] that B n is dominated by the quasi-real Compton singularity. In 
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Figure 14: B n for e'p — > e~p as function of the cm. scattering angle at different 
beam energies, as indicated on the figure. The calculations are for different hadronic 
intermediate states (X) in the blob of Fig. [12] : N (dashed curve), irN (dashed-dotted 
curves), sum of the N and ttN (solid curves). The data points are from the A4 Coll. 
(MAMI) [99] . Calculations and figure from Ref. [T08] . 



this (extreme forward limit) case, the hadronic tensor can be expressed in terms of the 
total photo-absorption cross section on the proton, aZ, allowing to express B n through 
the simple analytic expression : 



B r , 
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i Q2 

mi 



(61) 



One notices that the quasi-real Compton singularity gives rise to a (single) logarithmic 
enhancement factor which is at the orgin of the relatively large value of B n . 

In Fig. [151 the estimate from Ref. [HOj based on Eq. (|6T!) is shown for different 
parameterizations of the total photo-absorption cross section. The beam normal spin 
asymmetry has been measured at SLAC (E-158) at an energy E e = 46 GeV (-y/i ~ 
9 GeV) and very forward angle (Q 2 ~ 0.05GeV 2 ). First result |102j indicate a value 
B n ~ —3.5 — > —2.5 ppm, confirming the estimate shown in Fig. [151 

At intermediate energies, around E e ~ 3 GeV, and forward angles, B n has also been 
measured by the HAPPEX and GO Collaborations. The simple "diffractive" formula of 
Eq. (16"TI) does not rigorously apply any more and one has to calculate corrections due to 
the deviation from forward scattering. Such calculation have recently been performed 
in Refs. |110[ I112j in different model approaches, where the calculation of Ref. |112| 
includes sub leading terms in Q 2 . The predicted asymmetries are in basic agreement with 
first results reported by HAPPEX [Iffi] and GO [THS] . 

In Table [21 the present status of measurements of beam normal spin asymmetries, 
most of them in the forward angular range, is shown. 



4.3.4 Hard scattering region 

In the hard scattering region, the beam and target normal spin asymmetries were esti- 
mated in Refs. [1051 l3Tj through the scattering off a parton, which is embedded in the 
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Figure 15: Energy dependence of B n for e~*p —>■ e~p (right panel) at Q 2 = 0.05 GeV 2 
(corresponding with very forward scattering angle) using Eq. flSTl) for different param- 
eterizations of the total photo-absorption cross section on the proton (left panel) : 
Block-Halzen log fit (dotted curves), Block-Halzen log 2 fit (solid curves), and Donnachie- 
Landshoff fit (dashed curves). Calculations from Ref. [HOj . 

Table 2: Summary of measurements of the beam normal spin asymmetry in elastic 
electron-proton scattering 



EXP. 


E e (GeV) 


Q 2 (GeV 2 ) 


B n (ppm) 


SAMPLE [98] 


0.192 


0.10 


-16.4 ± 5.9 


A4 [99] 


0.570 


0.11 


-8.59 ± 0.89 


A4 [99] 


0.855 


0.23 


-8.52 ± 2.31 


HAPPEX [100] 


3.0 


0.11 


-6.7 ± 1.5 


GO |103| 


3.0 


0.15 


-4.04 ± 1.05 


GO |103| 


3.0 


0.25 


-4.81 ± 2.03 


E-158 |102j 


46 


0.06 


-3.5 -> -2.5 



nucleon through a GPD. 

The GPD estimate for the target spin asymmetry A n for the proton is shown in the 
left-hand plot of Fig. [16] as a function of the CM scattering angle for fixed incoming 
electron lab energy, taken here as 6 GeV. Also shown is a calculation of A n including the 
elastic intermediate state only [30]. The result, which is nearly the same for either of 
the two GPD parameterizations which were used in Ref. [37], is of order 1%. Fig. [T6l on 
the right also shows a similar plot of A n for a neutron target. The predicted asymmetry 
is of opposite sign, reflecting that the numerically largest term in Eq. ( |58i) is the one 
proportional to Gm- The results are again of order 1% in magnitude, though somewhat 
larger for the neutron than for the proton. A precision measurement of A n is planned at 
JLab [113j on a polarized 3 He target; it will provide access to the elastic electron-neutron 
single-spin asymmetry from two-photon exchange. 

Using the same phenomenological parametrizations for the GPDs, B n was found to 
yield values around +1 ppm to +1.5 ppm in the few GeV beam energy range, see Fig.[T71 
In particular, the forward angular range for e~^p — > e~p scattering was found to be a 
favorable region to get information on the inelastic part of B n . Because in the handbag 
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Figure 16: Nucleoli analyzing power A n , which is equal to the normal polarization P n . 
The elastic contribution is shown by the dotted curve [30] . The GPD calculation for the 
inelastic contribution is shown by the dashed curve for a gaussian GPD model, and by 
the solid curve for a modified Regge GPD parameterization. The GPD calculation is 
cut off in the backward direction at — u = M 2 . In the forward direction it goes down to 
Q 2 = 2 GeV 2 (modified Regge GPD) and to Q 2 = M 2 (Gaussian GPD). Figure from 



calculation, real and imaginary parts are linked, a direct measurement of B n may yield 
a valuable cross-check for the real part, which was found crucial in understanding the 
unpolarized cross section data for e~p — > e~p at large momentum transfer. 
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Figure 17: B n for elastic e~p scattering as function of e at different values of Q 2 as 
indicated on the figure. The upper thick curves (B n > 0) are the GPD calculations for 
the kinematical range where s, — u > M 2 . For comparison, the elastic contribution is also 
displayed : lower thin curves (B n < 0). Calculations and figure from Ref. |105j . 
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5 Conclusions and open issues 



The striking difference between the unpolarized (Rosenbluth) and polarization transfer 
measurements of the proton Ge/Gm form factor ratio has triggered a renewed interest in 
the field of two-photon exchange in electron-nucleon scattering experiments. Theoretical 
calculations both within a hadronic and partonic framework, which were reviewed in 
this work, made it very likely that hard two-photon exchange corrections are the main 
culprit in the difference between both experimental techniques. Despite the long history 
of two-photon exchange corrections, it is interesting to note that concepts developed 
over the past decade, such as generalized parton distributions which describe two-photon 
processes with one or two large photon virtualities, enter when quantifying two-photon 
exchange corrections at larger Q 2 . 

The mo del- independent finding is that the hard two-photon corrections hardly affect 
polarization transfer results, but they do correct the slope of the Rosenbluth plots at 
larger Q 2 in an important way, towards reconciling the Rosenbluth data with the Ge/Gm 
ratio, which decreases linearly in Q 2 , from the polarization data. This being said, it is 
fair to state at this point that neither of the current calculations convincingly quantifies 
the effect either due to uncontrolled approximations (such as using ad hoc assumptions 
for off-shell vertices in hadronic calculations), or extrapolations beyond their region of 
validity (such as applying partonic calculations in the low momentum transfer regime). 
The present and forthcoming high-precision electron scattering data, which aim at testing 
two-photon exchange effects, have shifted the emphasis, however, from the qualitative to 
the quantitative realm. 

Besides entering as corrections to elastic electron-nucleon scattering data, two-photon 
exchange processes are also the leading corrections to improve our quantitative under- 
standing of the hyperfine splitting (hfs) in hydrogen. We reviewed in this work that our 
best estimate, based on the recent data for polarized nucleon structure functions which en- 
ter the polarizability correction to the hydrogen hfs, yields a result A po i = (1.3±0.3) ppm. 
Combined with recent estimates of the Zemach radius, the correction to the hydrogen 
hfs falls short of the data by about 2 ppm, or about 2.8 standard deviations. 

We have also reviewed the recent measurements of large beam normal single spin 
asymmetries (SSA), of order —10 to —100 ppm, and measured with ppm level accuracy, 
which arose over the past few years as an interesting spin-off of the high precision parity 
violation experiments in electron-nucleon scattering. We discussed how such asymme- 
tries, which measure the imaginary part of the two-photon exchange amplitude, can be 
expressed through non-forward nucleon structure functions and provide a new tool to 
access hadron structure information. 

We end this review by spelling out a few open issues and challenges (both theoretical 
and experimental) in this field : 

1. Experimental measurements of the two-photon exchange processes 

In order to use electron scattering as a precision tool, it is clearly worthwhile to 
arrive at a quantitative understanding of two-photon exchange processes. This 
calls for detailed experimental studies, and several new experiments are already 
planned. A first type of experiments is to perform high precision unpolarized ex- 
periments and look for nonlinearities in the Rosenbluth plots. A recent study |114] 
found that the non-linearities in e in current Rosenbluth data stay small over a 
fairly large Q 2 range. Such experiments however cannot separate that part of the 
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two-photon corrections which is linear in e and which appears to be the dominant 
source of corrections in theoretical calculations. The difference between elastic e~ 
and e + scattering on a proton target directly accesses the real part of the two- 
photon exchange amplitude and the part linear in e. At intermediate Q 2 values 
some forthcoming e + /e~ experiments at JLab [50] and VEPP-3 |115] will allow one 
to systematically study and quantify two-photon exchange effects. Also forthcom- 
ing is an experimental check of the predicted small effect of two-photon processes 
on the polarization data by measuring the e dependence in polarization transfer 
experiments |49j . 

2. Further refinements on the theoretical side : real part of the two-photon exchange 
amplitude 

The calculations of the real part of the two-photon amplitude for elastic eN scat- 
tering described in this work are performed either within a hadronic framework at 
lower Q 2 values or partonic picture at larger Q 2 values. The hadronic calculation 
assumes some off-shell vertex functions and performs a four- dimensional integral 
of the box diagram. It would be desirable to have a dispersion relation framework 
where the real part is explicitely calculated as an integral over the imaginary part, 
which can be related to observables (i.e. on-shell amplitudes) using unitarity. In 
the forward limit such dispersion relations are underlying, e.g., the evaluation of 
the hyperfine splitting in hydrogen. In the non-forward case, the convergence of 
such dispersion relations would need further study. Within the partonic framework, 
a full calculation of the processes with two- active quarks (the so-called cat's ears 
diagrams) remains to be performed. 

3. Systematic exploration of the imaginary part of the two-photon exchange amplitude 

The ongoing experiments for both the beam and target normal SSA in elastic eN 
scattering will trigger further theoretical work, and a cross-fertilization between 
theory and experiment can be expected in this field. 

4. Two-photon exchange effects in resonance production processes 

The sensitivity of unpolarized measurements of Ge/Gm t° two-photon exchange 
effects also signals that such potential effects can be expected when measuring reso- 
nance transition form factors out to larger Q 2 values as is performed, e.g., at JLab. 
In particular a first theoretical study of the two-photon exchange contribution to 
the eN — > eA(1232) — > enN process with the aim of a precision study of the 
ratios of electric quadrupole (E2) and Coulomb quadrupole (C2) to the magnetic 
dipole (Ml) 7*iVA transitions, has been recently performed |116j . The two-photon 
exchange amplitude has been related to the iV — > A GPDs. It was found that 
the C2/M1 ratio at larger Q 2 values depends strongly on whether this quantity 
is obtained from an interference cross section or from the Rosenbluth-type cross 
sections, in similarity with the elastic, eN — > eN, process. It will be interesting to 
confront these results with upcoming new Rosenbluth separation data at interme- 
diate Q 2 values in order to arrive at a precision extraction of the large Q 2 behavior 
of the Rem and Rsm ratios. 

5. Two-photon exchange effects in deep-inelastic scattering processes 

Two-photon exchange effects were also studied in the past in deep-inelastic scat- 
tering processes; see, e.g., Refs. |117l 11181 1119] . They can be expected to affect 
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the extraction of the longitudinal structure function from L/T Rosenbluth type 
experiments in the deep-inelastic scattering region, which remains to be quantified. 

6. Two-photon exchange effects in parity violating elastic electron-nucleon scattering 
and related study of^yZ box diagram effects 

In recent years an unprecedented precision has been achieved in parity violating 
electron scattering experiments. Two-photon exchange processes may be relevant 
in the interpretation of the current generation of high precision parity-violation 
experiments, in particular for using those data to determine the strange-quark 
content of the proton; for a first study see Ref. [59] . 

A related issue in the interpretation of parity violating electron scattering exper- 
iments are the 7Z box diagram contributions. Such calculations have been per- 
formed for atomic parity violation [120] corresponding with zero momentum trans- 
fer as well as in the deep-inelastic scattering region [121] . by calculation the 7Z 
exchange between and electron and a quark. A full calculation in the small and 
intermediate Q 2 regime, where many current parity violation experiments are per- 
formed is definitely a worthwhile topic for further research. 

7. Two-photon exchange effects in hydrogen hyperfine splitting 

More accurate data, and data at lower Q 2 , is forthcoming on polarized inelastic 
e-p scattering. These data will allow a more precise evaluation of the polarizability 
corrections to hydrogen hfs. In addition, a measurement of hfs in muonic hydrogen 
may be possible [22] • The calculated corrections for muonic hydrogen have different 
weightings because of the muon mass, and with calculations such as [HDl EU ES] 
updated with new scattering data, the result of the measurement could be presented 
as an independent accurate determination of the proton's Zemach radius. 
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